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ABSTRACT 

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study 
a foreground quasar's environment in absorption. We used a sample of 17 Lyman limit systems with 
column density TVhi > 10^^ cm~^ selected from 149 projected quasar pair sightlines, to investigate the 
clustering pattern of optically thick absorbers around luminous quasars at z ~ 2.5. Specifically, we 
measured the quasar-absorber correlation function in the transverse direction, and found a comoving 
correlation length of — %.2t_\^ h^^ Mpc (comoving) assuming a power law correlation function, 
^ oc r~'^ , with 7 — 1.6. Applying this transverse clustering strength to the line-of-sight, would 
predict that 15 — 50% of all quasars should show a A'hi > 10^^ cm~^ absorber within a velocity 
window of Au < 3000 km s^^. This overpredicts the number of absorbers along the line-of-sight by 
a large factor, providing compelling evidence that the clustering pattern of optically thick absorbers 
around quasars is highly anisotropic. The most plausible explanation for the anisotropy is that the 
transverse direction is less likely to be illuminated by ionizing photons than the line-of-sight, and that 
absorbers along the line-of-sight are being photoevaporated. A simple model for the photoevaporation 
of absorbers subject to the ionizing flux of a quasar is presented, and it is shown that absorbers with 
volume densities n-n < 0.1 cm~'^ will be photoevaporated if they lie within ^ 1 Mpc (proper) of 
a luminous quasar. Using this simple model, we illustrate how comparisons of the transverse and 
line-of-sight clustering around quasars can ultimately be used to constrain the distribution of gas in 
optically thick absorption line systems. 

Subject headings: quasars: general - intergalactic medium - quasars: absorption lines - cosmology: 
general - surveys: observations 
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1. INTRODUCTION 

Although optically thick absorption line systems, that 
is Lyman Limit Systems (LLSs) and damped Lyman-a 
systems (DLAs), are detected as the strongest absorp- 
tion lines in quasar spectra, the two types of objects, 
quasars and absorbers, play rather different roles in the 
evolution of structure in the Universe. The hard ultravi- 
olet radiation emitted by luminous quasars gives rise to 
the ambient extragalactic ultraviolet (UV) background 
(see e.g. Haardt & Madau 1996; Miralda-Escude, 2002 
iMeiksinI 12005) responsible for maintaining the low neu- 
tral fraction of hydrogen 10~^) in the intergalactic 
medium (IGM), established during reionization. How- 
ever, high column density absorbers represent the rare lo- 
cations where the neutral fractions are much larger. Gas 
clouds with column densities log A^hi > 17.2 are optically 
thick to Lyman continuum (tll ^ 1) photons, giving rise 
to a neutral interior self-shielded from the extragalactic 
ionizing background. In particular, the damped Lya sys- 
tenis dominate the neutral gas content of the Universe 
(jProchaska et a l. 2005), which is the primary reservoir 
for the formation of stars observed in local galaxies. 

One might expect optically thick absorbers to be 
absent at small separations from luminous quasars. For 
example, a quasar at z = 2.5 with magnitude r — 19, 
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emits a flux of ionizing photons that is 130 times higher 
than that of the extragalactic UV background at an 
angular separation of 60", corresponding to a proper 
distance of 340 kpc. Indeed, the decrease in the 
number of optically thin absorption lines (log A'hi < 17.2 
hence tll ^ 1), in the vicinity of quasars, known 
as the proximity effect (Baithk et al. 1 988). provides 
a me asurement of the UV background ijScott et alJ 
|2000|) . If Nature provides a nearby background quasar 
sightline, one can also study the transverse proximity 
effect, which is the expected decrease in absorption 
in a background quasar's Lya forest, caused by the 
transverse ionizing flux of a foreground quasar. The 
transverse effect has yet to b e detected, in spite of 
many atte mpts llCrottsI Il989t iDobrzvcki fc Bechtolc 
1991 ; iFernand ez-Soto. Barcons . Carballo. fc Webb 

19951 ^iskc & Will iger, I2OO1F 

Schirber. Miralda-Escude. fc Mcnonaldl I200I fcTofll 



been 



2004 but see Jakobsen etaL 2003). 

On the other hand, it has long 
known that quasars are associated with en- 
hancements in the distribution of galaxies 



jBahcall, Schmidt. & Gunn' "196?; ' Yee fc Gre en Tgsl 
1987;_BahcaU fc Chokshi.l99L: .Smith Bovle. fc Maddoxl 
200firBrown. Bovle. fc WebsteHl200lUSerber et a,lJl2006t 
Coil et al.M200Q') . although these measurements of quasar 
galaxy clustering are mostly limited to low redshifts 
< 1.0. Recentlv. rAdelberger fc Steidell ((2005), measured 
the clustering of Lyman Break Galaxies (LBGs) around 
luminous quasars in the redshift range (2 < 2: < 3.5), and 
found a best fit correlation length of tq = 4.7 Mpc 
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Fig. 1. — Distribution of foreground quasar redshifts, trans- 
verse separations, and ionizing flux ratios probed by tlie back- 
ground quasar sightlines. The upper plot shows the ratio of the 
ionizing flux to the UV background guVi note the general 
trend. The lower plot shows foreground quasar redshift versus 
transverse comoving separations. The (blue) squares have a Keck 
spectrum of the background quasar, (red) triangles have Gemini 
background spectra, (magenta) upside down triangles have MMT 
background spectra, and (green) circles have SDSS background 
spectra. Filled symbols outlined in black have a super-LLS within 
|Ai;| < 1500 km s~^ at the foreground quasar redshift (see Ta- 
ble Q and open symbols have no absorber. The region to the left 
of the dotted line is excluded by the SDSS flber collision limit of 
6 = 55", which explains the paucity of SDSS background spectra 
there. The Keck/Gemini/MMT spectra probe angular separations 
an order of magnitude smaller than the flber collision limit, allow- 
ing us to study the foreground quasar environment down to down 
to 70 kpc where the ionizing flux is ~ 10,000 times the UV 
background. 



(7 = 1.6), very si milar to the auto-correl at ion length of 
z ^ 2 - 3 LBGs (|Adelberger et al.ll200^ . ICooke et al.l 
l)2006|) recently measured the clustering of LBGs around 
DLAs and measured a best fit tq = 2.9 Mpc 
with 7 = 1.6, but with large unce rtainties (sec also 
iGawiser et aD l200lt iBouche fc Lo^enthal .2004.) . If 
LBGs are clustered around quasars, and LBGs are 
clustered around DLAs, might we expect optically 
thick absorbers to be clustered around quasars? This is 
especially plausible in light of recent evidence that DLAs 
arise from a high redshift galaxy populatio n which are 
not unlike LBGs l|Schavell2n0lt iMoller et al.l l2002). 

Clues to the clustering of optically thick absorbers 
around quasars come from proximate DLAs, which have 
absorber redshifts within 3000 km s~^ of th e emission 
redshi ft of the quasars (see e.g. lM0ller et aliri99&). Re - 
cently. iRusseh et al.l()2006j) fsee also lEllison et al.ll20"0^ . 



compared the number density of proximate DLAs per 
unit redshift to th e average number dens ity of DLAs in 
the the Universe ijProchaska et aljl2005j) . They found 
that the abundance of DLAs is enhanced by a factor of 
~ 2 near quasars, which they attributed to the clustering 
of DLA-galaxies around quasars. 

Projected pairs of quasars with small angular separa- 
tions {9 <5') but different redshifts, can also be used to 
study the clustering of absorbers around quasars. In this 
context, clustering is manifest as an excess probability, 
above the cosmic average, of detecting an absorber in the 
neighboring background quasar spectrum, near the red- 
shift of the foreground quasar. In the first of this series 
of four papers on op tically thick absorbers near quasars 
ijHennawi et al.! '2006. henceforth QPQl), we used back- 
ground quasar sightlines to search for optically thick ab- 
sorption in the vicinity of foreground quasars: 149 pro- 
jected quasar pairs were systematically surveyed for Ly- 
man limit Systems (LLSs) and Damped Lya systems 
(DLAs) in the vicinity of 1.8 < z < 4.0 luminous fore- 
ground quasars. A sample of 27 new absorbers were un- 
covered with transverse separations R < 5 h^^ Mpc (co- 
moving) from the foreground quasars, of which 17 were 
super-LLSs with A^^hi > 10^^ cm'^. 

The distribution of foreground quasar redshifts, trans- 
verse separations, and ionizing flux ratios probed by the 
projected pair sightlines studied in QPQl are illustrated 
in Figure The filled symbols outlined in black indi- 
cate the sightlines which have an absorption line sys- 
tem with A^Hi > 10^^ cm~^ within a velocity interval of 
\Av\ = 1500 km s~^ of the foreground quasar redshift, 
and the open symbols represent sightlines with no such 
absorber. The line density of absorbers per unit redshift 
at z ^ 2.5 w ith column densities TVht > 10^^ cm^^ is 
dN/dz ~ O.S llO'Meara et all 1200^ . which imphes that 
the expected number of random quasar-absorber coinci- 
dences within |Ati| = 1500 km s""'^ is ^ 3%. The ex- 
pected number of random quasar-absorber associations 
in the sample of 149 sightlines is (N) — 4.4; whereas, 17 
systems were discovered in QPQl - compelling evidence 
for clustering. 

This clustering is particularly conspicuous on small 
scales: out of 8 sightlines with R < 500 kpc, 4 
quasar-absorber pairs were discovered, which implies a 
small scale covering factor of ~ 50%, whereas {N) = 0.18 
would have been expected at random. This excess of 
~ 25 over random may not come as a surprise when 
one considers that galaxies are strongly clustered around 
quasars on small scales. However, the light from every 
isolated quasar in the Universe traverses these same small 
scales on the way to Earth. Naively, we would then ex- 
pect a comparable covering factor of A^hi > 10^^ cm~^ 
super-LLSs along the line-of-sight (LOS) - which is def- 
initely not observed to be the case! 

How can we quantify the clustering of absorbers around 
quasars? How do we compare the transverse clustering 
pattern to the number of proximate absorbers observed 
along the LOS? What can the quasar-absorber cluster- 
ing pattern teach us about quasars and absorbers? These 
questions are addressed in this second paper of the se- 
ries. We briefly review the QPQl survey of projected 
quasar pairs and discuss the quasar-absorber sample in 
§ El A formalism for quantifying the line-of-sight (LOS) 
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and transverse clustering of absorbers around quasars 
is presented in § |3 In § 0] we introduce a maximum 
likelihood technique for estimating the quasar-absorber 
correlation function, which we use to measure the clus- 
tering of our quasar-absorber sample in § El Our mea- 
surement of the transverse clustering is used to predict 
the expected number of proximate absorption line sys- 
tems which should be observed along the line-of-sight 
in § El In § 13 we introduce a simple analytical model 
for the photoevaporation of optically thick absorbers by 
quasars, which is used to illustrate how comparisons of 
the line-of-sight and transverse clustering can constrain 
the distribution of gas in optically thick absorbers. We 
summarize and discuss our results in jj |H1 

Paper III of this series (jHennawi fc Prochaskall200^ 
investigates fluorescent Lya emission from our quasar- 
absorber pairs, and echelle spectra of several of the 
qu asar-LLS sys tems published here are analyzed in Paper 
IV l|Prochaska fc Hcnnawi 2006). Throughout this paper 
w e use the bes t fit WMAP (only) cosmological model 
of lSpereel et"al . (2003), with Qrn = 0.270, = 0.73, 
h = 0.72. Unless otherwise specified, all distances are 
comoving*. It is helpful to remember that in the chosen 
cosmology, at a redshift of z = 2.5, an angular separa- 
tion of Ad = 1" corresponds to a comoving transverse 
separation of i? = 20 h^^ kpc, and a velocity difference 
of 1500 km s^^ corresponds to a radial redshift space 
distance of s = 15 Mpc. For a quasar at z = 2.5, 
with an SDSS magnitude of r = 19, the flux of ioniz- 
ing photons is 130 times higher than the ambient ex- 
tragalactic UV background at an angular separation of 
60" (comoving R = 1.2 Mpc). Finally, we use the 
term optically thick absorbers and LLSs interchangeably, 
both referring to quasar absorption line systems with 
logiVni > 17.2, making them optically thick at the Ly- 
man limit (tll ^ !)■ 

2. QUASAR-ABSORBER SAMPLE 

In this section we briefly summarize the quasar- 
absorber sample and the parent sample of projected pairs 
from which it was selected; for full details see QPQl. In 
§ O we will quantify the clustering of absorbers around 
quasars in the transverse direction by comparing the 
number of quasar-absorber pairs discovered to the num- 
ber from random expectation. Although lower column 
density systems were published in QPQl, we focus here 
on absorbers with iVni > lO^^cm"^, or so called super- 
LLSs. This choice reflects a variety of concerns. First, we 
wanted to maximize the number of quasar-absorber pairs 
- had we restricted consideration to only DLAs we would 
have been left with 5 systems. Second, this is the lowest 
column density at which we believe we can identify a sta- 
tistical sample of absorbers with the moderate resolution 
spectra used in QPQl. Finally, the column density dis- 
tribution of optically thick absorption line systems has 
been determined for logA^ni > 19 ( Peroux et al. 2005; 
lO'Meara et alJ20T)6l) : whereas, it is unknown at the lower 
column densities 17.2 < logiVni ^19. 

Modern spectroscopic surveys select against close pairs 
of quasars because of fiber collisions. For the Sloan Dig- 
ital Sky Survey (SDSS), the finite size of optical fibers 
implies only one quasar in a pair with separation < 55" 

* Note that in QPQl proper distances were quoted. 



can be observed spectroscopically on a given plate^, and a 
slightly smaller limit (< 30") applies for the Two Degree 
Field Quasar Survey (2QZ) ijCroom et alJl200^ . Thus, 
for sub-arcminute separations, additional spectroscopy 
is required both to discover companions around quasars 
and to obtain spectra of sufficient quality to search for 
absorption line systems. For wider separations, projected 
quasar pairs can be found directly in the SDSS spectro- 
scopic quasar catalog. Hcnnawi et al. (2006a) used the 
3.5m telescope at Apache Point Observatory (APO) to 
spectroscopically confirm a large sample of photometri- 
cally selected close quasar pair candidates, and published 
the largest sample of projected quasar pairs in existence. 

In QPQl we combined high signal-to-n oise ratio (SNR) 
modera te resolution spectra of the closest Hcn nawi et alJ 
l)2006aj) projected pairs, obtained from Gemini, Keck, 
and the Multiple Mirror Telescope (MMT), with a large 
sample of wider separation pairs, from the SDSS spec- 
troscopic survey. The Keck, Gemini, and MMT spectra 
were all of sufficient signal-to-noise ratio to detect ab- 
sorbers with column densities A'hi > lO^^cm"^, and a 
SNR criterion was used to isolate the SDSS projected 
pair sightlines for which such an absorber could be de- 
tected in the background quasar spectrum. We consid- 
ered all projected quasar pair sightlines which had a co- 
moving transverse separation of R < 5 Mpc, a red- 
shift separation Aw > 2500 km s~^ between the fore- 
ground and background quasar (to exclude physical bi- 
naries) , and which satisfied our SNR criteria, and arrived 
at a total of 149 projected pair sightlines in the redshift 
range 1.8 < z < 4.0. 

A systematic search for optically thick absorbers with 
redshifts within |A?;| < 1500 km of the foreground 
quasars was conducted by visually inspecting the 149 
projected pair sightlines, where the velocity window was 
chosen to bracket the uncertainti es in the f oreground 
quasars' systemic redshift (see e.g. lRichards~et al. 200'^ 
and § 4 of QPQl). Voigt profiles were fit to systems with 
significant Lya absorption to determine the H I column 
densities. We uncovered 27 new quasar absorber pairs 
with column densities 17.2 < logA^ni < 20.9 and trans- 
verse comoving distances 71 h^^ kpc < R < 5 Mpc 
from the foreground quasars, of which 17 were super- 
LLSs with logiVni > 19. The mean redshift of the fore- 
ground quasar in the parent sample was (z) = 2.47 and 
that of the 17 quasar-super-LLS pairs was (zabs) = 2.55. 

The completeness and false positive rate of the 
QPQl survey a re a significant source of concern. 
iProchaska et alJ l)2005t) demonstrated that the spectral 
resolution (FWHM ~ 150 km s"^) and SNR of the 
SDSS spectra are well suited to constructing a complete 
(> 95%) sample of DLAs (logA^m > 20.3) at z > 2.2; 
however, the completeness at the lower column densities 
logA^Hi > 19, considered here has not been systemati- 
cally quantified. Furthermore, aggressive SNR criteria 
were employed in QPQl in order to gather a sufficient 
number of projected quasar pairs. Line-blending can sig- 
nificantly depress the continuum near the Lya profile and 

^ An exception to this rule exists for a fraction {~ 30%) of 
the area of the SDSS spectroscopic survey covered by overlapping 
plates. Because the same area of sky was observed spectroscopi- 
cally on more than one occasion, there is no fiber collision limita- 
tion. 
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TABLE 1 

SuPER-LLSs Near, Quasars from QPQl 



Name 


Zbg 


Zfg 


Ae R 

(") (ft-i kpc) 


Zabs 


|Ad| 
(km s^^) 


Avf^ 
(km s ■'-) 


log TVhi 
(cm-2) 


guv 


Telescope 


SDSSJ0225— 0739 


2.99 


2.440 


214.0 


4310 


2.4476 


690 


500 


19.55 ± 0.2 


5 


SDSS 




3.14 


2.308 


3.7 


72 


2.3025 


540 


1500 


20 45 + 2 


6369 


Keck 




3.55 


3.387 


179.0 


4195 


3.387 


20 


1000 


10 2^ + 025 


20 


SDSS 


SDSSJ0338-0005 


3.05 


2.239 


/ O.O 


1 /1 1 c; 
1410 


2.2290 


960 


1500 


20.9 ±0.2 


13 


SDSS 


SDSSJ0800+3542 


2.07 


1.983 


23.1 


415 


1.9828 


40 


300 


19.0 ±0.15 


488 


Keck 


SDSSJ0833+0813 


3.33 


2.516 


103.4 


2112 


2.505 


980 


1000 


19.45 ± 0.3 


18 


SDSS 


SDSSJ0852+2637 


3.32 


3.203 


170.9 


3917 


3.211 


550 


1500 


19.25 ±0.4 


13 


SDSS 


SDSSJ1134+3409 


3.14 


2.291 


209.2 


4073 


2.2879 


320 


500 


19.5 ±0.3 


11 


SDSS 


SDSSJ1152+4517 


2.38 


2.312 


113.4 


2216 


2.3158 


370 


500 


19.1 ±0.3 


30 


SDSS 


SDSSJ1204+0221 


2.53 


2.436 


13.3 


267 


2.4402 


370 


1500 


19.7 ±0.15 


625 


Gemini 


SDSSJ1213+1207 


3.48 


3.411 


137.8 


3246 


3.4105 


30 


1500 


19.25 ±0.3 


39 


SDSS 


SDSSJ1306+6158 


2.17 


2.111 


16.3 


302 


2.1084 


200 


300 


20.3 ±0.15 


420 


Keck 


SDSSJ1312+0002 


2.84 


2.671 


148.5 


3129 


2.6688 


200 


500 


20.3 ±0.3 


23 


SDSS 


SDSSJ1426+5002 


2.32 


2.239 


235.6 


4529 


2.2247 


1330 


500 


20.0 ±0.15 


19 


SDSS 


SDSSJ1430-0120 


3.25 


3.102 


200.0 


4517 


3.115 


960 


1500 


20.5 ±0.2 


26 


SDSS 


SDSSJ1545+5112 


2.45 


2.240 


97.6 


1873 


2.243 


320 


500 


19.45 ±0.3 


30 


SDSS 


SDSSJ1635+3013 


2.94 


2.493 


91.4 


1861 


2.5025 


820 


500 


> 19 


111 


SDSS 



Note. — Optically thick absorption line systems near foreground quasars. The background and foreground quasar 
redshifts are denoted by z^^^ and Zfg, respectively. The angular separation of the quasar pair sightlines is denoted by 
A6, which corresponds to a transverse comoving separation of R at the foreground quasar redshift. Absorber redshift is 
indicated by z^bsi Eind \Av\ is the velocity difference between the absorber redshift and our best estimate of the redshift 
of the foreground quasar. Our estimated error on the foreground quasar redshift is denoted by Aiifg. Foreground quasar 
redshifts and redshift errors were estimated according to the detailed procedure described in § 4 of QPQl. The logarithm 
of the column density of the absorber from a fit to the H I profile is denoted by log A'^hi ■ The column labeled "Telescope" 
indicates the instrument used to observe the background quasar. The quantity guv = 1 + ^Qso/^UVB is the maximum 
enhancement of the quasars ionizing photon flux over that of the extragalactic ionizing background at the location of the 
background quasar sightline, assuming that the quasar emission is isotropic (see Appendix A of QPQl). We compare to 
the UV background computed by F. Haardt & P. Madau (2006, in preparation) 



mimic a damping wing, biasing column density measure- 
ments high or giving rise to false positives. Any statis- 
tical study will thus suffer from a 'Malmquist'-type bias 
because line-blending biases lower column densities up- 
ward, and the line density of absorbers dN / dz, is a steep 
function of column density limit. 

Based on visually inspecting the 149 background 
quasar spectra and a comparison with echelle data for 
three systems, we estimated that the QPQl survey 
was ~ 90% complete for logiVni > 19.3 for all the 
Keck/Gemini/MMT spectra and ~ 3/4 of the SDSS 
spectra, which accounts for about 125 of the 149 spectra 
searched. To address the false positive rate, we com- 
pared with echelle data for three systems and found that 
our column density was overestimated by ~ 2.5ct for one 
system, raising it above the super-LLS (logA^Hi > 19) 
threshold. However, this absorber was located blueward 
of the quasars Ly/3 emission line, in a part of the spec- 
trum 'crowded' by the presence of both the Lya and Ly/3 
forests. A more careful examination of the completeness 
and false positive rate of super-LLSs identified in spec- 
tra of the resolution and SNR used in QPQl is definitely 
warranted. In § we explore how our clustering mea- 
surement changes if we discard systems within w Icr of 
the logA'^Hi = 19 threshold. 

Relevant quantities for the super-LLS-quasar pairs 
which are used in our clustering analysis are given in 
Tabled The distribution of foreground quasar redshifts, 
transverse separations, and ionizing flux ratios probed by 
all of our projected pair sightlines is illustrated in Fig- 
ure Q The filled symbols outlined in black indicate the 



sightlines which have an absorption line system near the 
foreground quasar with logiVni > 19 (see Table ^| and 
open symbols are sightlines with no such absorber. 

3. QUANTIFYING QUASAR-ABSORBER CLUSTERING 

In the absence of clustering, the line density of ab- 
sorption line systems per unit redshift above the column 
density threshold A^hi is given by the cosmic average 

/ dN\ c 
(^)(>A^Hi,.) = nA/e..^, (1) 

where n is comoving number density of the galaxies or 
objects which give rise to absorption line systems, A is 
their absorption cross section (in comoving units), fcov 
is the covering factor, and H{z) is the Hubble constant. 
Note that the line density is degenerate with respect to 
the combination nAfcov and only their product can be 
determined by measuring the abundance of absorption 
line systems. 

At an average location in the Universe, the probability 
of finding a absorber in a background quasar spectrum 
within the redshift interval Az = 2(1 + z)Av/c, corre- 
sponding to a velocity interval 2Av is simply P = {dN / 
dz)Az. For a projected pair of quasars, clustering around 
the foreground quasar will increase the probability of 
finding an absorber in the background quasar spectrum 
in the vicinity of the foreground quasar; whereas, the 
foreground quasars radiation field could reduce this en- 
hancement by photoevaporating absorption systems. If 
the quasar sightlines have a comoving transverse separa- 
tion R, and assuming that one searches a velocity inter- 
val ±Au about the foreground quasar redshift (because 
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1 10 100 

Tabs (proper kpc/h) 

Fig. 2. — Dependence of line-of-sight correlation function on the 
(proper) size of the absorber cross section Tabs- A velocity interval 
of Av = 3000 km s~^ was assumed. Small cross-sections sample 
smaller scales near the quasar where the correlation ^qa oc 
function is large, giving a larger value for the average in eqn. J3 
X||; whereas, large cross-sections dilute ^qa in the average, giving 
smaller values of The solid (black) line is for a correlation 

length of r = 4.7 Mpc and slope of 7 = 1.6. The dotted (red) 
line is for a correlation function with slope 7 = 2.0, normalized to 
the have the same clustering amplitude as the the shallower 7 = 1.6 
model at r = 1 Mpc. 

of redshift uncertainties) , we can express the increase in 
line density near the foreground grounds in terms of a 
transverse correlation function x±{R) a-s 

dN /dN\ , . , 

az \ az / 

where x± {Rj ^1^) is given by an average of the 3-d quasar- 
absorber correlation function, S,qa{t), over a cylinder 
with volume V — A[2Av/aH{z)]. Here, a is the scale 
factor and 2Av{l + z)/H{z) is the length of the cylinder 
in the LOS direction. We can thus write 

X±{R)^Ip l^dV ^QA{r) (3) 

where Z is a comoving distance in redshift space. The 
last approximation in eqn. ^ assumes that the volume 
average over the cylinder can be replaced by the line 
average along the line-of-sight direction, which is valid 
provided that we are in the 'far-field' limit, i.e., the trans- 
verse separation is much larger than the diameter of the 
cylinder R ^ \/A. Thus provided we consider distances 
R much larger than the dimension of the absorber, the 
transverse clustering is independent of the absorption 
cross-section. 

For proximate absorbers along the line of sight, we can 
similarly define a line-of-sight correlation function 

dN/dz = (dN/dz) [1 + x\\ (Aw)] (4) 

where 

X\\{Av) = f^dZ f dA^QA(r). (5) 



The lower limit of the radial integration is set to Zcut, a 
cutoff which is introduced to parametrize our ignorance 
of the geometry of the absorbers. For instance, if the 
absorbers were pancake shapes which were always ori- 
ented perpendicular to the line of sight (like face-on spi- 
ral galaxies) then we would not cutoff the line-of-sight 
integration at all {Zcut = 0), since face on pancakes at 
zero separation can still obscure the quasar. For 'hard 
sphere', we would set Zcut = 2-y/A/7r — presuming that 
the obscuring cross-section of the absorber drops to zero 
for points interior to it. As smaller values of Zcut will 
correspond to larger line-of-sight clustering, we hence- 
forth conservatively assume the hard sphere case and use 

Zcut = 2^ Aj-K 

Note that one is no longer in the 'far-field' limit for 
the integral in eqn. 0, and the clustering amplitude X|| 
explicitly depends on the cross section of the absorber. 
It is easy to understand the nature of this dependence if 
one considers that the line density cx nA is fixed by mea- 
surements of {dN/dz). Thus, smaller cross-sections cor- 
respond to larger volume number densities, but smaller 
scales near the quasar are being sampled by the integral 
in eqn. jSJl, where the correlation ^qa cx; r~'^ function 
is large. Conversely, larger cross-sections sample regions 
further from the quasar, and thus ^qa is averaged over 
a larger volume and hence diluted. Figure |21 illustrates 
the dependence of X||(Aw) on the proper radius of the 
absorber cross section, Tabs, for a range of sizes. 

In eqns. Q and Q CQA(f ) represents a real space cor- 
relation function and it may appear that we have ne- 
glected redshift space distortions. Strictly speaking, the 
redshift space correlation function ^qa{R, Z), should be 
appear in eqns. ^ and [3 This £,qa{R, Z) is the convolu- 
tion of the real space correlation function ^(r), with the 
velocity distribution in the radial direction, which can 
have contributions from both peculiar velocities and un- 
certainties in the systemic redshifts of the quasar. How- 
ever, provided that the distance in redshift space over 
which we project Aw contains most of the probability 
under this distribution, it is a good approximation to re- 
place the redshift space correlation function, under the 
integrals in eqns. I@J and ©, with the real space correla- 
tion function, because radial velocities will simply move 
pairs of points within the volume. 

4. ESTIMATING THE CORRELATION FUNCTION 

4.1. Maximum Likelihood Estimator 

Given a quasar-absorber correlation function i^qa, 
eqns. |(5J) and I^J describe how to compute the prob- 
ability of finding an absorber in the redshift interval 
Az = 2(1 -I- z)Av/c at a transverse distance R from a 
foreground quasar: P{R, z) = {dN/dz) Az. Considering 
that we only have 17 quasar-absorber pairs selected from 
149 sightlines, it will be difficult to measure more than 
a single parameter with reasonable errors. Hence, we as- 
sume the quasar-absorber correlation function to have a 
power-law form 

eQA = c(^) \ (6) 

where C is a clustering amplitude and rg is the cor- 
relation length. The amplitude C is degenerate with 
To, but we choose to estimate C, because it allows 
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for the possibility of anti-correlation (C < 0) which 
could result from the QSO ionizing radiation field. Mo- 
tivated by the th e slope of the L BG auto-correlation 
function (Adclbcr ger et al.l l2005a(l . we choose to fix 
7 = 1.6. A similar pr ocedure was employed by 
l|Adelberger fc Steidell \20()H. henceforth AS05), who 
measured clustering of LBGs around luminous quasars 
(2 ^ ^ ^ 3.5), and found a best fit correlation length of 
To = 4.7 Mpc. We set tq = 4.7 Mpc as a fiducial 
value; thus C can be interpreted as the quasar-absorber 
clustering amplitude relative to the AS05 quasar-LBG 
result. 

Consider an ensemble of N projected pair sightlines 
with background quasars at transverse separations Ri 
from foreground quasars at redshifts Zi. Given ^qa, we 
can compute the associated probabilities Pi = P{Ri,Zi). 
Suppose that A/'slls of these A^ sightlines show absorp- 
tion from super-LLSs with logA^Hi > 19. The likelihood 
of the data, given the model parameter C is then 

WSLLS W-A^SLLS 

^{c)= n n (1-^^)' 

i 3 

where the probabilities Pi <1 are capped at 1. By max- 
imizing the likelihood with respect to the parameter C, 
we estimate the clustering amplitude from the data. 

4.2. Monte- Carlo Simulations 

We must verify that the maximum likelihood estimator 
in eqn. lO is unbiased, and we would also like to know 
how to assign error bars to an estimate of C. Both of 
these points can be addressed with monte-carlo meth- 
ods. The distribution of redshifts and transverse sepa- 
rations in Figure n is not uniform, and it is thus impor- 
tant that we preserve this distribution when constructing 
mock data sets to assign errors. 

For a 'true' value of the clustering amplitude C, we can 
compute the probabilities Pi = P{Ri,Zi) of observing 
an absorber for each of the N projected pair sightlines. 
Mock data sets can then be constructed by generating 
an A^-dimensional vector of deviates from the uniform 
distribution Xi, and assigning sightlines with Xi < Pi an 
absorption line system. The same maximum likelihood 
estimator is applied to these mock data sets to deter- 
mine an estimate C for each one, thus allowing us to 
measure the probability distribution P{C\C) of the esti- 
mate about the true model. 

4.3. Binned Correlation Function 

An alternative to the maximum likelihood technique 
described above, would be to estimate the correlation 
length by fitting a power-law model to estimates of 
the correlation function in bins of transverse separation. 
However, because we have only 17 quasar-absorber pairs, 
the best-fit correlation length from this procedure would 
be very sensitive to the binning chosen. We compute the 
correlation function in bins because it provides a useful 
way to visualize the data, although we do not fit the 
binned data for C 

A simple estimator for the transverse correlation func- 



tion in eqn. (01 is ^ 

where (QA) is the number of quasar-absorber pairs in a 
transverse radial bin centered on R, and (QR) — J2i {'^-^ I 
dz)Az, is the number of quasar-random pairs expected. 

4.4. The Column Density Distribution 

Before we can estimate the correlation amplitude with 
eqn. Q, we require the cosmic average line density of 
super-LLSs {dN/dz) . The line density is the zeroth mo- 
ment of the column density distribution 




[>N^i,z)= / fm{Nm.z)dNm- (9) 



lO'Meara et al.l f2006') measured the column density dis- 
tribution for super-LLSs in the range 19 < logA^Hi < 
20.3 at z ~ 2.7 and found good agreement with a power 
law f{N) = A{Nm/lQ^^ cm'^)'', where A = 1.10x10-^° 
and b = -1.43. For logAfm > 20.3, prochaska et alJ 
I2005j) measured the column density distribution from the 
SDSS quasar sample in redshift bins. 

To evaluate the i ntegral in eqn. we use the 

lO'Meara et al.l l|200fiD power law fit in the range 19 < 
log A^H T < 20.3, and a spline fit to the IProchask a et alJ 
l)2005j) results in the redshift bin centered at z ~ 2.7. 
For the redshift evolution, we simply scale the line den- 
sity of all absorbe rs by the evolution of D LAs log A^hi > 
20.3 measured bv IProchaska et aTl 1)20051) (see their Fig- 
ure 8). This procedure assumes that the abundance of 
super-LLSs evolves similarly to that of DLAs, an as- 
sumption that is con sistent with but not confirmed by 
lO'Meara et"all l|200^ . 

5. CLUSTERING RESULTS 

We applied the maximum likelihood estimator to the 
17 quasar- absorber pairs (see Tabled]) selected from the 
149 projected pair sightlines shown in Figure Q The 
maximum likelihood value of the clustering amplitude is 
C = 2.9, which corresponds to a correlation length of 
ro = 9.2 /i-i Mpc. 

In Figure |31 we show the probability distribution of 
the clustering amplitude maximum likelihood estimates 
P(C), for a model with true value equal to our mea- 
surement of C = 2.9. This distribution was created by 
applying the maximum likelihood estimator in eqn. Q 
to 100,000 mock realizations of our 149 projected pair 
sightlines, as described in § 14.21 The mean of this dis- 
tribution is (C) = 2.94, the dispersion is cr = 0.86, and 
the (red) dotted curve shows a Gaussian with the same 
mean and dispersion. The monte-carlo simulation indi- 
cates that our estimator is indeed unbiased and that the 
error distribution is reasonably well approximated by a 
Gaussian. Using the full distribution from the monte- 
carlo simulation, we compute the 68% confidence inter- 
val about the 'true' input value, which we use to quote 
errors on our clustering measurement: C — 2.9 ± 0.8, 

^ This is actually an estimator of x± averaged over the 

volume of the bin. We ignore this distinction, because the dif- 
ferences are substantially smaller than our errors and we are not 
fitting parameters from this estimate. 
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Fig. 3. — Distribution of clustering amplitude estimates from 
monte-carlo simulations. The black histogram shows the proba- 
bility distribution of the clustering amplitude maximum likelihood 
estimates P(C), for a model with true value equal to our measure- 
ment of C = 2.9. This distribution was created by applying the 
maximum likelihood estimator in eqn. J7J to 100,000 moc k rea liza- 
tions of our 149 projected pair sightlines, as described in i; l4.2l The 
mean of this distribution is (C) = 2.94, the dispersion is cr = 0.86, 
and the (red) dotted curve shows a Gaussian with the same mean 
and dispersion. The full distribution (black histogram) is used to 
assign errors to our measurement. 



T ' ' ' — ' — ' — ' ' ' I 




R (kpc/h) 

Fig. 4. — Binned transverse quasar-absorber correlation function 
for 17 quasar-absorber pairs selected from 149 projected pair sight- 
lines. The lines indicate the best-fit model from our maximum like- 
lihood analysis and the shaded regions indicate the range allowed 
by la errors estimated from monte-carlo simulations. The solid line 
and the blue shaded region are for an assumed correlation function 
slope of 7 = 1.6. The steeper dotted line and red shaded region are 
for 7 = 2.0. The dashed curve and green region show the transverse 
correlation function x± if we set the quasar-absorber correlation 
function ^QA to the QSO-LBG correlation function measured by 
AS05 (ro = 4.7 ± 1.3 /i-l Mpc; 7 = 1.6) 



or ro = 9.2t|^ y h ^ Mpc. If we assume a steeper 
slope for the correlation function 7 = 2, we obtain 
ro = 5.8+J:g h-^ Mpc. 

We detect strong clustering of super-LLSs around 
quasars. Figure ?? shows the binned transverse corre- 
lation function x± compared to the two maximum likeli- 
hood fits (7 = 1.6 and 7 = 2) and to the prediction from 
the quasar-LBG correlation function measured by AS05. 



For 7 = 1.6, the quasar-absorber clustering amplitude is 
three times larger than the quasar-LBG measurement of 
AS05 and inconsistent with it at the > la level. 

The auto-correlation functions of high redshift galax- 
ies tend to become progressively steeper on (comoving) 
scales R < 1 Mpc chara.cteristic of the sizes of 
dark matter ha los ijCoil et all 120051 lOuchi et alJ 120051 
iLee et aP l2006t ?), thus similar behavior might be in 
cross-co rrelation around quasars. Indeed, Hcnnawi et alJ 
l)2006a|) detected an order of magnitude excess quasar 
auto-clustering on scales < 100 h^^ kpc. For this rea- 
son, we also quote our results for 7 = 2, although a 
steeper slope is not clearly favored by the binned data in 
Figure^ 

At first glance, it may seem odd that the relative error 
on our measurement of ro is ~ 20% using only 17 quasar- 
absorber pairs; whereas AS05 quote a ~ 30% error on ro 
from ~ 200 quasar-LBG pairs around a sample of ^ 50 
quasars. The estimator used here is qualitatively simi lar 
to that used by AS05 (see also'Adelberecr et al. 2005a) in 
that both techniques rely on the radial clustering of pairs 
of objects at fixed angular positions, because the angula r 
selection functions are unknown (see IAdelbergerll2005bt) . 
It is thus worth explaining that our smaller errors arise 
from a few factors. 

First, the relative error decreases with clustering am- 
plitude, and our best fit has a factor of three larger 
clustering strength. Second, AS05 exclude scales R < 
1.2 Mpc from their analysis; whereas four out of our 
17 quasar-absorber pairs have R < 1.2 Mpc. Small 
scale pairs are effectively 'worth' many large scale pairs 
because the signal to noise per pair is much higher. Fi- 
nally, the AS05 errors are determined from the field to 
field dispersion in the data, which includes a contribu- 
tion from cosmic variance. Because the projected pairs 
of quasars are distributed over the entire sky, our mea- 
surement does not suffer from cosmic variance errors. 

To illustrate that the errors from our monte-carlo 
technique are sensible and comparable to the AS05 re- 
sult, we randomly resampled our pair sightlines with 
R > 1.2 Mpc to create a mock projected pair sam- 
ple 30 times larger, but with the same distribution of 
redshifts and transverse distances. We increased the 
search window Av — 3000 km s^^, in closer agree- 
ment with the I = 30 Mpc radial window averaged 
over by AS05. The average number of quasar- absorber 
pairs expected from this hypothetical enlarged sample 
is (N) = J2i Pi = 206, which is of order the number 
of LBG-AGN pairs used by AS05. Assuming 7 = 1.6 
and ro = 4.7 Mpc, our monte-carlo simulation gives 
ro = 4.7ti;o Mpc, or a relative error of ~ 20%, com- 
parable to the AS05 errors but slightly smaller. 

5.1. Systematic Errors 
5.1.1. 'Malmquist' Bias 

As discussed in § |2 the identification of logA'^Hi = 19 
in spectra at the SNR and resolution used in QPQl, 
can result in a significant 'Malmquist' type bias 
because line-blending scatters lower column density 
absorbers upward, and the line density of absorbers 
dN/dz, is a steep function of column density limit. 
If absorbers with column densities logiVni < 19 
scatter up into our sample, this error would bias 
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our clustering measurement high. To investigate the 
impact of this bias, we redo the clustering analysis 
but ignoring the six quasar-absorber pairs in Table ^ 
which have column densities within la of the threshold 
logNm = 19: SDSSJ 0256-^0039, SDSSJ 0800-f3542, 
SbSSJ 0852-1-2637', SDSSJ 1152-f4517, 

SDSSJ 1213-^1207, and SDSSJ 1635-f3013. The 
eleven remaining quasar-absorber pairs give a maximum 
likelihood clustering amplitude of C = 1.7 ± 0.7 or 
ro ~ 6.4j]![ g h^^ Mpc, compared to our measurement 

of C = 2.9 ± 0.8, or tq = 9.2j:J j h-^ Mpc for the 
full sample. This illustrates that even under very 
conservative assumptions about the possible effect of 
Malmquist bias, the clustering amplitude would be only 
be reduced by ~ 1.5(7. 

5.1.2. Redshift Errors 

Another possible source of bias in our sample could 
arise from the determination of the foreground quasar 
redshifts. Note that the clustering analysis properly 
takes the redshift uncertainties into account, by aver- 
aging the correlation function over a window of |Aw| = 
1500 km s~ . However, assigning redshifts to the fore- 
ground quasar (see § 4 of QPQl) is not a completely 
objective process. The quasar emission lines sometimes 
exhibit mild-BAL or metal line absorption, and the line 
centering can be sensitive to how these features are 
masked. It is possible that we were biased towards in- 
cluding quasar-absorber pairs in our sample, and thus 
tended to assign redshifts resulting in velocity differ- 
ences |Af I < 1500 km s~^. To address this issue we 
redo the analysis discarding the three quasar-absorber 
pairs in Table H which have velocity differences larger 
than the quoted error, Avfg > \Av\: SDSSJ 0225-0739, 
SDSSJ 1426-^5002, and SDSSJ 1635+3013. The 14 
remaining quasar-absorber pairs give a maximum like- 
lihood clustering amplitude of C = 2.4^Q y or ro = 

8.1^7^ Mpc. Furthermore, if we discard the five 
quasar-absorber pairs with the largest velocity differ- 
ences, we measure C = 1.9±0.7 or rg = 7.1^17 Mpc. 
Thus we conclude that even if we were biased towards 
including several quasar-absorber pairs with the velocity 
differences |At;| > 1500 km s~^, this would only reduce 
the the measured clustering amplitude by ~ la. 

6. ANISOTROPIC CLUSTERING OF ABSORBERS AROUND 

QUASARS 

If the clustering pattern of optically thick absorbers 
around quasars is isotropic, then we can use our estimate 
of the quasar- absorber correlation function in §010 pre- 
dict the number of 'proximate' super-LLSs that should 
be observed in a given velocity window. Aw, along the 
line-of-sight. According to eqn. the clustering en- 
hancement x\\ can be computed by averaging ^qa over a 
cylinder of length Av/aH{z) and cross sectional area A. 

However, because the cross-section of the absorbers is 
unknown, relating the transverse clustering to the line-of- 
sight clustering requires an assumption about the size of 
the cross section. Very little is known about the sizes of 
DLAs and LLSs. iBriggs et all |i989) detected H I 21cm 
absorption of a DLA with A'hi = 5 x 10^^ cm^^ against 
an extended radio source, allowing them to place a lower 
limit of r > 12 h^^ kpc on the size of the absorbing 



region. IProchaskal l)1999(l estimated an absorbing path 
length of ~ 2h~^ kpc for a super-LLS with logTVni = 
19.1, by comparing the column density to the total vol- 
ume density, which was determi ned from the coU isionally 
excited C II* A1335 transition. lAdelbergerlTalJ (|200fll 
constrained the size of a DLA with logiVni = 20.4 to 
he > 8 kpc, based on their detection of fluores- 
cent Lyo! emission from the edge of a DLA-galaxy sit- 
uated ^ 8 kpc from the location of DLA absorption. 
iLopez et al.l ll2005j) measured identical column densities 
of TVhi = 102°-5 cm-2 for both DLAs at (z = 0.9313) de- 
tected in the individual images of the gravitational lens 
HE 0512 — 3329, allowing them to const rain the size of 
the DLA to be > 4/1"! kpc lSmette et all (see also^995). 
It is not clear how to relate these measurements to the 
absorption cross sections of the logiVni > 19 of interest 
to us here. 

To determine the cross section siz e as a functio n of 
hmiting column, rabs(> A^Hi) = \/ aA{> TVhO/tt, we 
adopt the simple approximation that the comoving num- 
ber density of absorption line systems, n, and the cov- 
ering factor, /cov are independent of column density, 
which gives the simple scaling Tabs oc a/ {dN/ dz) from 
eqn. Two families of sizes, 'small' and 'large' are 

considered, which we believe bracket the range of pos- 
sibilities. Fiducial physical sizes of Tabs = 5 kpc 
and Tabs — 20 h^^ kpc are chosen at the DLA thresh- 
old (logTVHi > 20.3) for the 'small' and 'large' ab- 
sorbers, respectively. The Tabs oc y/ {dN/ dz) scaling 
then predicts respective sizes of Tabs = 9 h^^ kpc and 
Tabs = 38 kpc for super-LLSs with logA'^Hi > 19. 

In Figure we show the probability that a quasar at 
z = 2.5 will have a 'proximate' optically thick absorption 
line systems within Av < 3000 km s~^, as a function of 
limiting column density. The left (right) panel shows the 
prediction for the 'small' ('large') family of cross-section 
sizes, and the dot-dashed curves shows the prediction in 
the absence of clustering (x|| = 0), which is simply the 
integral of the column density distribution in eqn. 
Since we have measured the transverse clustering only 
for logA^Hi > 19, this figure assumes that clustering is 
independent of limiting column density. 

The transverse clustering overpredicts the fraction of 
quasars which have a proximate absorber by a large fac- 
tor. For example, for 'small' absorption cross sections 
and 7 = 1.6, the quasar-absorber correlation function 
measured from the transverse direction predicts that a 
fraction P = 0.30 ± 0.07 of all quasars should show 
a proximate super-LLS (logA^ni > 19) within Av — 
3000 km s^ . Even the AS05 clustering amplitude, 
which is a factor of ~ 3 smaller than our best fit, would 
predict P — 0.12lg Qg. A steeper correlation function re- 
sults in even more proximate absorbers. For the best fit 
transverse clustering amplitude with 7 = 2, the proba- 
bility would be P = 0A9to\l. 

Making the absorbers larger changes this prediction 
by factors of ~ 2 — 3. The 'large' absorbers predict a 
fraction P = 0.14±0.03of quasars should have a proxi- 
mate absorber, given our best-fit clustering amplitude for 
7 = 1.6. Our steeper 7 = 2 fit gives P = 0.14^°;°^ and 
the AS05 result predicts P — 0.07 ± 0.03 and our steeper 
7 = 2 fit gives P = 0.14j;o °4. Although the line density 
of proximate super-LLSs near quasar has yet to be mea- 
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Fig. 5. — Probability of a quasar at 2 ~ 2.5 having a proximate absorber within 3000 km s^^ as a function of limiting column density. 
The left panel is for 'small' absorber cross sections and the right panel is for 'large' absorbers. Our measurement of the quasar-absorber 
correlation function ^qa from the transverse direction is used to predict the number of proximate absorbers along the line-of-sight. The 
curves and shaded regions show the predictions and la errors for the quasar-absorber correlation function ^qa estimated from the transverse 
direction in §|21 as well as the clustering strength measured by AS05. The linestyles and colors correspond to the same models as in Figure|4] 
The lower dotted line indicates the prediction in the absence of clustering, which is just the integral of the column density distribution in 
eqn.|5] The point with the error bar represents the lRusscll et al. 1 2006,') measurement for proximate DLAs (log Afjji > 20.3) at 2 ^ 3.4. 



sured (but see lProchter. Hennawi. fc Prochaskal200(iD . it 

is incontrovertible that 15 — 50% percent of quasars do 
not show a super-LLS within Av < 3000 km s~^ along 
t he line of sight . 

iRussell et al.l ()2006 f) recently measured the line den- 
sity of proximate DLAs (logA^ni > 20.3) with Aw < 
3000 km s"^ from 731 SDSS quasars. The median red- 
shift of their quasar sample was z = 3.335 and they 
found (dN/dz) = OMt^ H- The point with the error 
bar in Figure [3 indicates the probability of having a 
proximate DLA from the Russell et al. (2006) measure- 
ment, Pdla = 0.013^g;oo2- Our best-fit clustering am- 
plitude predicts that P = 0.12 ± 0.3 (P = 0.25j:|^;^^) of 
quasars should have a nearby DLA for small absorbers 
and 7 = 1.6 (7 = 2). Large absorbers change this pre- 
diction to P = 0.055 ± 0.013 (P = 0.069t[!:[J^i^) ^. 

The caveat should be included that our prediction for 
the number of proximate absorbers from the transverse 
clustering, is very sensitive to the small scale behavior 
of the correlation function. Specifically, using eqn. |5l 
to predict the line-of-sight clustering implicitly assumes 
that the power law model of the correlation function 
^ oc , is valid down the lower limit of the integra- 
tion, Zcut, which we have set to be the diameter of the 
absorbers. For super-LLSs, we used (proper) radii of 
fahs = 9 kpc and Tabs = 38 kpc, for small and 
large absorbers, respectively, corresponding to (comov- 
ing) Zcut — 63 and Zcut = 266 at z — 2.5. The small- 
est transverse separation of our projected pair sample is 
R — 72 h^^ kpc, and there are just 5 sightlines with R < 
300 h^^ kpc, three of which have an absorber (see Fig- 
ure^and Table^l. The transverse measurement assumes 

Comparing the IRussell et al.l 12003) measurement with our 
prediction ignores evolution in the clustering and the line den- 
sity between the mean redshift of our sample {z = 2.5) and theirs 
(z = 3.335). The line density changes by ~ 50% over this redshift 
range fProchaska et al. 2005), but this is small compared to the 
large discrepancy with prediction of the transverse clustering. 



a single power law from R = 70 h^^ kpc — 5 h^^ Mpc. 
One should bear in mind that only a handful of projected 
pairs probe the small scales {R < 300 h^^ kpc) that the 
line-of-sight clustering is very sensitive too; but, it is re- 
assuring that closest bin in Figure 01 is consistent with 
the power law fit. 

We have shown that the transverse clustering which we 
quantified in §|Sloverpredicts the abundance of proximate 
absorbers along the line-of-sight by a large factor ^ 4 — 
20, under reasonable assumptions about the sizes of the 
cross sections of these absorbers. The clustering pattern 
of absorbers around quasars is thus highly anisotropic. 
The most plausible explanation for this anisotropy is that 
the transverse direction is less likely to be illuminated 
by ionizing photons than the line-of-sight, and that the 
optically thick absorbers along the line-of-sight are being 
photoevaporated. We will discuss the physical effects 
which could give rise to this anisotropy in §|H1 Next, we 
introduce a simple model that provides physical insight 
into the problem of optically thick absorbers subject to 
the intense ionizing flux of a nearby quasar. 

7. PHOTOEVAPORATION OF OPTICALLY THICK 
CLOUDS 

The problem of an optically thick absorption line sys- 
tem exposed to the ionizing flux of a nearby luminous 
quasar is analogous to that of a neutral interstellar cloud 
being exposed to the ionizing radiation o f an OB star - a 
proble m which was first investigated bv lOort fc SpitzeJ 
i|1955D . iBertoldil l|1989t) classified the behavior of pho- 
toevaporating clouds based on their initial column den- 
sity and the ionization parameter at the location of the 
cloud, and developed an analytical solution to follow the 
radiation-drive n implosion pha se of the cloud. Interest- 
ingly, although iBertoldH ('1989') was primarily concerned 
with the fate of interstellar clouds near OB stars in H II 
regions, he commented briefly on the applicability of the 
same formalism to Lya clouds exposed to the ionizing 
flux of a quasar. 
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Fig. 6. — Regions of parameter space where an optically thick absorber can survive near a luminous quasar set by egn. IIUI The dashed 
(red), black (solid) and dotted (blue) curves correspond to r = 17, 19, & 21 mag, respectively. Optically thick absorbers above a curve can 
survive, whereas those below are photoevaporated. Left: Survival regions in the volume density-distance plane, assuming a neutral column 
of log A^HI = 20.3. Right: Survival regions in the neutral column density-distance plane, assuming a volume density of nj] = 0.1 cm~'^. 



7.1. Cloud Zapping 

Following 'Bertoldjl l)1989(l . we model an optically thick 
absorber as a homogeneous spherical neutral gas cloud 
with total number density of hydrogen hh , which is em- 
bedded in photoionized intergalactic medium at temper- 
ature T = 20, 000 K, corresponding to an isothermal 

sound speed q — 16.5 km s , where is the 

temperature in units of 20, 000 K. If the cloud is at a 
distance r from a luminous quasar which is emitting S 
ionizing photons per second, then the the ionizing flux at 
the cloud surface, Fi — S/Anr'^, will drive an ionization 
front into the neutral gas. If the flux i s large enough 
to make the conditions R-type fsee iSnitze r (1978.) for a 
discussion of how ionization fronts are classified) , an R- 
type ionization front will propagate through the cloud 
without dynamically perturbing the neutral gas. As the 
front propagates into the cloud, it ionizes an increas- 
ing column of neutral gas, steadily reducing the ionizing 
flux at the front, until the conditions become M-type, at 
which point the front will stall and drive a shock into the 
neutral upstream gas. This shock will implode the cloud 
and compress it until conditions become D-type, allow- 
ing the front to continue prop agating and establishing a 
steady photoevaporation flow ijBertol di 1989). 

The cloud will be 'zapped', or completely photoevapo- 
rated, if the ionizing flux is large enough relative to the 
column density, such that the entire cloud can be ionized 
in a recombination time. In this case, the R-type front 
will completely cross the cloud without stalling. Part 
of the cloud will remain neutral and a shock will form 
provided that 



4947V, 



-1 

H,20.3 



r 



1.1 X lO""* Tso^ 



< 1 



(10) 



where 7Vh,20.3 is the total hydrogen column in units of 
2Q20.3 Qjj^s^ g^j^^^ ^jjg ionization parameter is defined by 
F = Fi/uQC with 



F = 2.58 X 10"^ 5*56 r 



Mpc 



(11) 



where 5*56 is the ionizing flux in units of 10^^ s ^, rMpc 
is the physical distance in units of Mpc, and rifj,! is the 
total hydrogen number density in units of 1 cm^"^. The 
condition that the ionization front be R-type at the sur- 
face of the cloud is F > 1.1 x 10-'*T2o^^ 

Consider our fiducial example of a foreground quasar 
with r = 19 at an angular separation of A9 = 1' from an 
absorber corresponding to a transverse proper distance of 
485 kpc. The ionizing flux is enhanced by guv — 130 over 
the UV background and = 5.2. For a DLA with total 
hydrogen column of A^h — lO'^" '^ cm~^ at this distance, 
the cloud will survive {S < 1) provided nn > 0.27 cm~^, 
which would give F < 0.002. 

The left panel of FigureElshows the lower limits on the 
volume density of a DLA with neutral hydrogen column 
log A^Hi = 20.3, set by the condition for cloud survival, as 
a function of physical distance from a quasar at z = 2.5. 
The dashed (red), black (sohd) and dotted (blue) curves 
correspond to r = 17, 19, & 21 mag, respectively. In 
the region below each curve, the volume densities are 
too small and the clouds are photoevaporated; whereas, 
above the curves the clouds can survive. The right panel 
shows the lower limit on the neutral hydrogen column 
density as a function of distance, assuming a volume den- 
sity riH = 0.1. If DLAs with log A'hi = 20.3 have sizes in 
the range Tabs ~ 1 — 5 kpc their corresponding number 
densities are nn ^ A^H/j^abs — 0.01 — 1. Thus according 
to Fignre^optically thick absorbers with rin ^0.1 will be 
photoevaporated if they lie within ^ 1 Mpc of a luminous 
quasar. 

Note that the condition in eqn. (O considers the total 
column density log A'h, not the neutral column, log A^hi- 
Although DLAs are expected to be predominantly neu- 
tral A'hi ~ A^H, the ionic absorption lines typically ob- 
served in LLSs suggest that they are the photoionized 
analogs of DLAs (Prochaska 1999) and that their ioniza- 
tion state is determined by ionization equilibrium with 
the UV background. In eqn. (^U)), we have assumed 
that optically thick absorbers are spherical top-hat den- 
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Fig. 7. — Toy model prediction for the probability of a quasar at 
z ~ 3.4 with i = 19.5 having a proximate IDLA within 3000 km 
as a function of volume density. The upper x-axis shows self- 
shielding radius, which is the minimum distance from the quasar 
at which a DLA can survive for a given volume density. The 
curves and shaded regions show the predictions and Icr errors for 
the quasar-absorber correlation function estimated from the 
transverse direction in §|S| as well as the clustering strength mea- 
sured by AS05. The linestyles and colors correspond to the same 
models as in Figure |^ The long-dashed horizontal lines indicate 
the measurement and ztltr range measured bv lRus sell et al. (2006) 
^'dla = 0-013^g'QQ2 . The short dotted horizontal line is the 
cosmic average (i.e. no clustering) measured by IProchaska et all 
120051) . This model suggests that DLAs with volume densities in 
the range 10~^ cm""^ < njj ^ 10~^ cm ~^ are required to a gree 
with the proximate DLA measurement of lR.ussell et al.i 120031 . 



sity distributions. Hence, an LLSs with column density 
TVhi ^ 20, can be thought of as a sightline that passes 
through the photoionized outskirts of a top-hat cloud 
which has a larger total hydrogen column A^h, and that 
the total column determines the ability of the cloud to 
survive a blast of ionizing radiation from the quasar. 
Thus for A^Hi ^ 20, we compute an ionization correc- 
tion with the standard approach, which is to assume a 
slab geometry and determine the ionization balance in a 
uniform ba ckgroun d, using a photoionization code such 
as Cloudy l|Ferland et al. 1998). To construct the curve 
in the right panel of Figure El wc interpolated through a 
grid of Cloudy solutions ^ for TVhiIA^h)- 

7.2. A Toy Model to Predict the Statistics of Proximate 

DLAs 

In this section we use a toy model to illustrate how 
quasar-absorber clustering can be used to constrain the 
physical properties of optically thick absorbers. The cri- 
teria i5 < 1 in eqn. H10|l . gives a minimum distance from 
the quasar, i?ss("-H) as a function of volume density, at 
which an optically thick absorption line system with col- 
umn density A^hi can survive. Our approach is to simply 
assume that absorbers at smaller distances are photoe- 
vaporated and absorbers at larger distances survive. We 
also assume that the transverse direction is not illumi- 

* The photoionization models used Cloudy version 6.0.2 and used 
a Haardt & Madau (2006, in preparation) QSO-|-galaxies spectrum 
at z = 2.5 and assumed a metallicity of -1.5 and J912 = —21.2. 



nated by the quasar, and hence the transverse cluster- 
ing measures the intrinsic quasar-absorber clustering, in 
the absence of ionization effects. Because proximate ab- 
sorbers are definitely illuminated, this intrinsic clustering 
is then reduced along the line-of-sight by photoevapora- 
tion. 

For the column density range logA^ni > 20.3, we eval- 
uate Rss at the lower limit logA^m = 20.3, which is a 
decent approximation because the column density distri- 
bution is steep, and the statistics will be dominated by 
absorbers near the threshold. To simplify the computa- 
tion, we take Rss to be a distance only along the line 
of sight, which is valid provided that Rss ^ VA. Then 
wc can write that the line density of absorbers within 
Aw < 3000 km s"^ is 

dN /dN\\ aRssHjz) , 1 

where x\\ is given by eqn. ||SJ) but with Z^ut — Rss- The 
assumption that the transverse direction gives the intrin- 
sic clustering in the absence of ionization effects amounts 
to using the correlation function ^qa, measured from 
the transverse clustering, in the line-of-sight integral in 
eqn. (Q. 

In Figure we show our toy model prediction for the 
probability of a quasar having a proximate DLA with 
log Ahi > 20.3 within Av < 3000 km s~^, as a function 
of total volume density of hydrogen rtn- We assumed 
that the DLAs have a size of Tabs — 5 kpc charac- 
teristic of the 'small' absorbers discussed in § |H| How- 
ever, identical results are obtained for 'large' absorbers 
(''abs = 20 kpc). Because our toy model excludes 
small scales Tabs ^ -Rss from the clustering integral in 
eqn. (jSJ, we are again in the far- field limit where the 
the volume average is nearly independent of the absorber 
cross section. This independence breaks down for very 
large densities tzh ^ 10^ j where i?ss approaches the size 
of the absorbers (see Figure [7|). The curves and shaded 
regions show the predictions and la errors for the trans- 
verse clustering measured in § as well as the clustering 
strength measured by AS05. We assumed a quasar with 
magnitude i — 19.5 at z = 3.335, chosen to match the 
mean magnitude and redsh ift of the proximate DLA sam- 
ple of IRussell et al.l lj2006D . The long-dashed horizontal 
lin es indicate the measu rement and ±lcr range measured 
bv lR.ussell etHI ll200fil) . Pdla = 0.013l^;gEJ^ . 

Our toy model suggests that DLAs with volume den- 
sities in the range 10^^ cm^'^ ^ n^i < 1 0~^ cm~'^ are 
required to agree with the measurement of 'Russell et alJ 
l|2006). Even the weaker AS05 clustering of LBGs around 
quasars would require nn ^ 1 cm~^. For larger volume 
densities, DLAs illuminated by quasars would survive at 
smaller radii where the clustering is strong, giving rise to 
a larger number of proximate DLAs. N ote that the com- 
parison of our toy model to the lRusseU et al. (2006) mea- 
surement in Figure 13 assumes that the clustering of ab- 
sorbers is independent of column density threshold (i.e. 
our measurement is for logA^ni > 19) and that there is 
negligible evolution in the clustering between the mean 
redshift of our clustering sample, z — 2.5 and that of the 
iR.ussell et a l. (2006) sample z — 3.335. These assump- 
tions were made because our quasar-absorber sample did 
not have sufficient statistics to measure the clustering 
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fsee ^[5]) in the sam e redshift and column density range as 
the lRussell et al.l (|^06) proximate DLA measurement. 

Although crude, this model illustrates how a com- 
parison of line-of-sight and transverse quasar ab- 
sorber clustering can be used to determine the den- 
sity distribution in optically thick absorbers. De- 
tailed models of self-shielding with radiative transfer 
llZheng k Miralda-Escud^l2002at iCantalupo et alJl200a 
IKollmeier et al.ll2006j) would be required for a more ac- 
curate treatment, and such analyses could easily include 
cuspy density profiles or a multiphase distribution of 
gas. Although better statistics and more theoretical work 
are necessary, the clustering of optically thick absorbers 
around quasars will provide important new constraints 
on the physical nature of these systems. 

8. SUMMARY AND DISCUSSION 

8.1. Summary 

In this paper we used a sample of 17 super-LLSs 
(log A^Hi < 19) selected from 149 projected quasar pairs 
sightlines in QPQl to investigate the clustering pattern 
of optically thick absorbers around quasars. Based on 
this data, this paper presented the following results: 

1. A simple formalism was presented for quantifying 
the clustering of absorbers around quasars in both 
the transverse and line-of-sight directions. The 
clustering of absorbers around quasars in the trans- 
verse direction is independent of the size of the 
absorption cross section; whereas, the line-of-sight 
clustering was shown to be sensitive to the cross 
section size. 

2. Applying this formalism to the 17 super-LLSs 
(logA^Hi < 19) selected from 149 projected quasar 
pair sightlines with mean redshift z = 2.5, we de- 
termined a comoving correlation length of tq = 
9.2^i['7 Mpc for a power law correlation func- 
tion with 7 = 1.6. This is three times stronger 
than the clustering of LBGs around quasars re- 
cently measured by AS05. If we assume a steeper 
slope of 7 = 2.0, we measure = S.S^J g Mpc. 

3. The clustering of optically thick absorbers around 
quasars is highly anisotropic. If we apply the clus- 
tering amplitude measured in the transverse direc- 
tion to the line-of-sight, the fraction of quasars 
which have a proximate absorber within A?; < 
3000 km s"^ is overpredicted by a factor as large 
as ~ 4 — 20, depending on assumptions about 
cross section sizes and the slope of the correlation 
function. The most plausible explanation for the 
anisotropy is that the transverse direction is less 
likely to be illuminated by ionizing photons than 
the line-of-sight, and that the optically thick ab- 
sorbers along the line-of-sight are being photoevap- 
orated. 

4. A simple model of absorbers as uniform spherical 
overdensities is discussed and we write down an 
analytic criterion which determines whether an ab- 
sorber illuminated by a quasar will be able to self- 
shield. This criterion indicates that optically thick 
absorbers with rifj ^0.1 will be photoevaporated 



if they lie within ^ 1 Mpc of a luminous quasar. 
We combine this criterion with a toy model of how 
photoevaporation affects the line-of-sight cluster- 
ing, to illustrate how comparisons of the line-of- 
sight and transverse clustering around quasars can 
ultimately be used to constrain the distribution of 
gas in optically thick absorption line systems. 

8.2. Discussion 

The anisotropic clustering pattern of absorbers 
around quasars suggests that the transverse direction 
is less likely to be illuminated by ionizing photons 
than the line-of-sight. This suggestion gains credi- 
bility in light of the null detections of the transverse 
proximity effect in t h e Lyg forests of projected 
quasa r p airs \C,mt\A rrnSflt inobrzvcki Rechtoldl 
|1991t iFernandez-Soto. Barcons. Carballo. fc WebH 
Il99a ILiske fc Willigeil pOOlF 

Schi rber. Miralda-Escude. fc McDonaldl l200l fCrofa 
l2004l but see Jakobsen et al.2003). Although these stud- 
ies are each based only on a handful of projected pairs, 
they all come to similar conclusions: the amount of (op- 
tically thin) Lya forest absorption, in the background 
quasar sightline near the redshift of the foreground 
quasar, is larger than average rather than smaller - 
the opposite of what is expected from the transverse 
proximity effect. Two physical effects can explain both 
the optically thin results and our result for optically 
thick systems: anisotropic emission or variability, which 
we discuss in turn. 

If quasar emission is highly anisotropic, the line-of- 
sight would be exposed to the ionizing flux of the quasar; 
whereas, transverse absorbers would be more likely to 
lie in shadowed regions. Studies of Type II quasars 
and the X-ray background suggest that quasars with lu- 
minosities comparable to our foreground quasar sample 

S Mr < —23) hav e ^ 30% of the solid angle obscured 
Ueda et all I200I iBarger et all l2f)0l iTreister fc TlTrl 
12005ft . although these estimates are highly uncertain. 
Naively, we would expect the covering factor of trans- 
verse absorbers to be approximately equal to the aver- 
age fraction of the solid angle obscured. But in QPQl we 
found a very high covering factor (6/8) for having an op- 
tically thick absorber with logTYni > 17.2 (see Figure 1 of 
QPQl) on the smallest (proper) scales R < 150 kpc. 
Although the statistics are clearly very poor, this high 
covering factor is suggestive of a significantly larger ob- 
scured fraction. 

If the ionizing flux of the foreground quasar varies con- 
siderably on a timescale shorter than the transverse light 
crossing time between the foreground and background 
sightlines, a transverse proximity effect might not be ob- 
servable. This is because the ionization state of the gas 
along the transverse sightlines is sensitive to the fore- 
ground quasars luminosity a light crossing time before the 
light that we observe was emitted. At 60" (1.2 h^^ Mpc) 
from a, z — 2.5 quasar the transverse light crossing time 
is 1.1 X 10^ yr. 

Currently, the lower limit on the intermittency of 
quasar emission comes from observations of the (op tically 
thin) proximity effect (Baitlik et al. 1988; Scott et all 
l200nj) . in the Lya forests near quasars. The presence of 
an optically thin proximity effect implies that the IGM 
has had time to reach ionization equilibrium with the 
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quasars increased ionizing flux, which requires that the 
duration of a burst of quasar radiation is lon ger than 
the I GM equihbration time, tburst ^ 10^ yr l)Martin3 
|2004() . The photoevaporation timescale for an optically 
thick absorber is ~ N-m/F or the light crossing time, 
whichever is longer, where F is the ionizing flux. At a 
distance of 100 kpc from an i = 19 quasar, it would take 
15, 000 yr to photoevaporate a DLA with logiVni = 20.3, 
i.e. comparable to the light crossing time if DLAs have 
sizes of '--^ 5 kpc. Thus the time it takes for the line-of- 
sight to manifest the effects of the quasars ionizing flux is 
~ 10"* yr for both the optically thin and optically thick 
regimes. Hence, if quasars emit in bursts of duration 
10^ yr ^ ^burst ^ 10^ yr, this would sufSce to explain the 
absence of the optically thin transverse proximity effect 
as well as the anisotropic clustering pattern of optically 
thick absorbers discovered here. 

The optically thin proximity effects are probably the 
most promising way to disentangle whether obscura- 
tion or variability can explain why the transverse di- 
rection is less likely to be illuminated by ionizing 
photons ISchirber. Miralda-Escude. fc Mc Donald 2004; 
lCroftl2004t) . This is because for the lower cohimn density 
absorbers characteristic of the Lya forest, cosmological 
N-body simulations can predict the statistical properties 
of the absorbers ah initio without the complications of 
radiative transfer. 

It might seem puzzling that we measure a clustering 
amplitude which is a factor of three larger than the clus- 
tering of LBGs around quasars at similar redshift mea- 
sured by AS05. One could argue that on small scales 
we are probing material intrinsic to the quasar. If in- 
terpreted as clustering, this intrinsic contribution would 
lead us to overestimate the correlation function. How- 
ever, only our closest sightline {R = 22 kpc, proper) 
is sufficiently small for this to be a worry, and this 
would imply a very large cross section for DLA absorp- 
tion. AS05 excluded small scales from their analysis, 
whereas, the strength of our signal is in part driven by 
the small scale systems. If we exclude sightlines with 
R < 500 kpc we measure a ro = T.Gtj'o'*"^ Mpc 
which is within la of the AS05 measurement. It is also 
possible that we measured a larger clustering signal be- 
cause the galaxies which host optically thick absorbers 
are more strongly biased with respect to the dark matter 
distribution than LBGs (but see Cooke ct al. 2006). 

We argued that the transverse clustering predicts that 
galaxies correlated with a quasar should give rise to 
an absorption line systems with logA'ni > 19 within 
Av < 3000 km s~^ a significant fraction of the time 
'--^ 15 — 50%. We postulated that these systems are not 
observed because of photoevaporation, but shouldn't a 
similar argument also apply to galaxies? High redshift 
galaxy spectra should show optically thick absorption 
due to nearby correlated galaxies, in addition to intrinsic 
absorption due to H I in the galaxy itself. The individual 
spectra of high redshift galaxies rarely have high enough 
signal-to-noise ratio to detect even DLAs. However, the 
break due to Lyman limit absorption can be detected in 
very deep integrations with the added advantage t hat the 
LLSs are much more abundant. Recently, Shaplcv et al.l 
l)2006(l detected emission blueward of the Lyman limit in 
2 out of 14 LBGs at z ~ 3. For a galaxy-absorber corre- 



lation length in the range ro = 5 — 10ft. ^ Mpc, eqn. ^ 
would predict that 25—75% of galaxies should show a cor- 
related absorber within a window of Av < 500 km s~^, 
comparable to the redshift uncertainties of LBGS. This 
number assumes LLSs with logiVui > 17.2 have a size of 
~ 50 kpc and we used th e abundance of LLSs mea- 
sured bv iPeroux et al.l l)2003j) . While it is possible that 
the lack of flux blueward of the Lyman limit in LBGs 
is due to absorption by gas intrinsic to the LBG, corre- 
lated absorbers along the line of sight could also play a 
signiflcant role. 

Using the ionizing flux of a quasar to study the dis- 
tribution of neutral gas in optically thick absorbers is 
a powerful new way to study these absorption line sys- 
tems. We showed how comparisons of the line-of-sight 
and transverse clustering of absorbers around quasars 
constrains their distribution of gas. Theoretical mod- 
els of LLSs and DLA s which in clude radiative trans- 
fer and self-shielding iZheng <^ Miralda-Escude 2002^ 
iCantalupo et al.ll2005t |Kollmeier et al. 2006) should ex- 
plore how morphology, cuspy density profiles, or a multi- 
phase medium, change the ability of absorbers to survive 
near quasars. 

Optically thick H I clouds in ionization equilib- 
riimr with a radiation field, re-emit 60% of the 
ionizing photons they absor b as fluorescent L y a re- 
combination line photons llGould fc Weinberg! 1996; 
Zheng fc Mirald a-Escude) '2002b t iCantalup o et al.l |200E : 
Adelberger et a l. 2006; KoUmeie ret al.l l2006'). Recently, 
Adelberger et a l. (2006) reported a detection of fluores- 
cence from a serendipitously discovered DLA (log iVni = 
20.4) situated 49" away from a luminous quasar ai z — 
2.84. Fluorescent Lya emission from absorbers illumi- 
nated by quasars offers a new window to study physi- 
cal properties of absorbers, as well as a second trans- 
verse sightline to observe the ionizing flux of a quasar. 
In QPQl we published seven new quasar-absorber pairs 
which would have fluorescent surface brightnesses of 
MLya = 19.5 — 24.3 mag arcsec"^, if the foreground 
quasars emit isotropically. However, the clustering 
anisotropy discussed here suggests that these transverse 
absorbers may not be illuminated, and the lack of a de- 
tection of fluorescent emission would provide compelling 
evidence in favor of this conclusion. Constraints on flu- 
orescent emission from these systems will be discussed 
in th e next paper in this series (Hcnnawi & Prochasl^ 
120061) . It is particularly intriguing that proximate DLAs 
observed along the line-of-sight, seem to preferentially 
exhibit Lya emission superimposed on the Lya ab sorp- 
tion trough l|M(aller et alJ 119981: lElIison et alJl2002il . Is 
this Lya emission fluorescence? Is the fluorescent emis- 
sion from proximate DLAs observable? 

In QPQl we argued that a factor of ^ 20 more trans- 
verse quasar-absorber pairs could be compiled in a mod- 
est amount of observing time, which would improve our 
measurement of the transv erse correlation func tion by a 
factor of ~ 5. The study of R.ussell et al.l lj2006i) searched 
~ 700 SDSS quasars for DLAs and found 12 proximate 
DLAs with Av < 3000 km s~\ However, the fuU 
SDSS quasar survey has ^ 10,000 quasars with z > 2.2 
((Schneider et al. 2006), which could be used to search 
for proximate systems and increase the number known to 
> 100. This would be sufficient to measure the column 
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density distribution /(A^athi) of absorbers near quasars, 
and any differences between it and the distribution in 
the 'field' would provide an important new constraint 
on the physical nature of optically thick absorbers. Fi- 
nally, we note that similar studies using projected quasar 
pairs and searching for line-of-sight proximate absorbers 
can be conducted for Mg II, C IV or other metal ab- 
sorption lin es systems near quasars fBowen ' et~aLll2?M 
IProchter. H cnnawi. & Prochaska 2006), yielding similar 
insights into their physical nature. 

Large samples of optically thick absorption line sys- 
tems near quasars are well within reach, for transverse 
systems using quasar pairs as well as proximate absorbers 
along the line-of-sight. This data will provide new oppor- 
tunities to characterize the environments of quasars, the 
physical nature of absorption line systems, and it will un- 
cover new laboratories for studying fluorescent emission 
from optically thick absorbers. 
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